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ABSTRACT 

A parametric model is presented for the blended- 
wing-body airplane, one concept being proposed 
for the next generation of large subsonic trans- 
ports. The model is defined in terms of a small 
set of parameters which facilitates analysis and 
optimization during the conceptual design process. 
The model is generated from a preliminary CAD 
geometry. From this geometry, airfoil cross sections 
are cut at selected locations and fitted with analytic 
curves. The airfoils are then used as boundaries for 
surfaces defined as the solution of partial differential 
equations. Both the airfoil curves and the surfaces 
are generated with free parameters selected to give 
a good representation of the original geometry. The 
original surface is compared with the parametric 
model, and solutions of the Euler equations for 
compressible flow are computed for both geometries. 
The parametric model is a good approximation of 
the CAD model and the computed solutions are 
qualitatively similar. An optimal NURBS approxi- 
mation is constructed and can be used by a CAD 
model for further refinement or modification of the 
original geometry. 
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NOMENCLATURE 

D 0 ,Di Dirichlet boundary conditions 
N 0 ,Ni Neumann boundary conditions 
P Point on NURBS surface 

P NURBS control points 

X Surface coordinate 

C Airfoil chord length 

F\ , F 2 Airfoil thickness parameter 

Mi , M 2 Airfoil camber parameters 

N NURBS basis function 

Pi , P 2 Airfoil camber parameters 

T Airfoil thickness parameter 

h Wing section span 

it, v Surface parametric variables 

w NURBS weight 

x,y,z Surface coordinates 

x, y Airfoil coordinates 

x e , y e Airfoil trailing edge 

y CAD airfoil coordinate 

y c Airfoil camber 

y t Airfoil thickness 

0 Twist angle 

£, r] Wing section parametric variables 

a PDE weighting factor 

1. INTRODUCTION 


Many of the busiest airports in Asia will run out of 
space in the near future. One solution to the problem 
proposed by the aircraft industry is the development 
of huge airliners carrying up to 800 passengers. 
Along with the tradition aircraft configurations, 
some manufacturers, including McDonald Douglas 1 , 
have advanced a flying wing configuration similar to 
the B2 bomber. This unconventional blended-wing- 
body (BWB) concept for future subsonic transports 
is an aerodynamically efficient airplane with the 
interior passenger volume distributed over a large 
portion of the centerbody. 
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New concepts in aircraft designs often begin with 
only a rough sketch on the designers note pad. 
These conceptual designs always have to be modified 
due to structural and aerodynamic constraints. 
For this reason, it is convenient to have a model 
which is defined by a set of parameter values so the 
the model can be easily modified as more design 
variables are considered. This report will describe 
the application of a conceptual design tool that will 
facilitate aerodynamic analysis in the early stages 
in the development of an aircraft. The parametric 
model which is generated in this manner now can 
be optimized to improve the aerodynamic charac- 
teristics of the aircraft by varying a small set of 
parameters, or the complete parameter space can be 
examined to determine all possible feasible designs. 
By performing a more complete analysis at the early 
stages of the development, the developer will be 
less likely to have to make major design changes in 
the final Computer-Aided Design (CAD) models to 
be used for Computational Fluid Dynamics (CFD) 
analysis and wind tunnel testing. 


ure 2. Four of these cuts were on the main horizon- 
tal section of the wing and two were on the vertical 
winglet. These cross sections were output as a dense 
set of points. The airfoils in the parametric model 
were then generated with the parameters selected to 
give the best fit to the CAD airfoils. The paramet- 
ric airfoil equations defining the coordinates x(£) and 
y(£) of the points on an airfoil are 


x (0 = Csimrt, y(£) = y t (0 + y c (0, 


Vt(0 


VciO = { 


T 

— — (sin2TT^ + FisinA-K^ + F 2 sin67r£), 
j£(2Pir — r 2 ) if r < Pi 

_^ F ( r _P 2 )2( 2r + p 2 _ 3 p 1 ) 

-(7^ F (r-Pi) 2 (2r + P 1 -3P 2 ) 
if Pi < r < P 2 

? T%r(l - m + 2 P 2 r - r 2 ) if r > P 2 , 


r = siring 

0 < Pi < P 2 < 1, 0 < f < 1. (1) 


This report will describe an application of the 
conceptual design tool Rapid Airplane Parametric 
Input Design (RAPID) 2,3 in creating a parametric 
model for the the BWB airplane. Since RAPID 
was developed for modeling conventional airplane 
geometries, modifications were necessary to apply 
the techniques to the current application. The 
input parameters for a geometry may come from 
a variety of sources and range from crude sketches 
to detailed drawings or surfaces generated from a 
CAD system. For the BWB geometry a preliminary 
CAD model was available and was used to gener- 
ate most of the parameters used in the RAPID model. 

2. AIRFOIL APPROXIMATION 

The construction of the parametric model for the 
BWB begins with the definition of airfoil cross sec- 
tions. The actual CAD data was used to define the 
airfoil parameters in RAPID. The number of airfoils 
defining the model determines the total number of 
parameters in the model and also effects the fidelity 
of the model. If a large number of cross sections 
are made, the CAD model can be more accurately 
approximated by the RAPID model. However, the 
number of parameters increases with the number of 
airfoils, and there is the ever present danger of os- 
cillations in the model that is inherent in this as 
well as many other interpolation and approximation 
schemes. From the CAD surface given in Figure 1, 
six cuts were made at the positions indicated in Fig- 


The design parameters in these equations allow suf- 
ficient freedom to give a good approximation to all 
the airfoil cross sections used for the model. The pa- 
rameter C is the section chord length. The thickness 
curve yt(0 is defined with the thickness parameter 
T and Fourier coefficients F\ and F 2 . A general- 
ization of the piecewise quadratic camber curve in 
the RAPID system has been developed to model the 
wing sections of the BWB configuration. The cam- 
ber curve y c (0 is a piecewise cubic curve which can 
have both a local maximum and a local minimum. 
The camber curve has relative extrema of M\ and 
M 2 at locations P\ and P 2 , respectively. The cam- 
ber locations are the fractions of total chord length 
measured from the airfoil trailing edge. If P± = P 2 
and Mi = M 2 this expression reverts to the original 
piecewise quadratic camber curve in RAPID. The ad- 
ditional parameters are needed for the BWB model 
to approximate the forward camber and aft reflex of 
the centerbody airfoil. 

The selection of a set of design parameters to fit a 
given airfoil can be accomplished by solving a non- 
linear optimization problem. Assume that a set of 
values i = 1 , • • • , m have been selected based on 
some desired criteria for distributing points around 
the airfoil. From these values, a corresponding set of 
values Xi, i = can be calculated from the 

first equation in (1). Now use the airfoils from the 
CAD data and calculate corresponding y coordinates 
yi , i = 1 , • • • , m, by interpolating the data if neces- 
sary. Observe that on the upper and lower surfaces 
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of the airfoil the value of y can be given as a single- 
valued function of x. The chord length C can be 
measured and will not be included in the set of vari- 
ables for the optimization problem. Thus, the best 
fit to the given CAD data is the airfoil defined by 
equation ( 1 ) with the remaining design parameters 
T, Fi , F 2 , Mi , M 2 , Pi , P 2 chosen to minimize 

m 

X>Cr, f 1 ,f 2 ,m 1 ,m 2 , p u p 2 ) - yif 

The solution of this nonlinear least-squares problem 
defines the airfoil shapes for the parametric model. 
This optimization problem can be easily solved by 
a number of methods. The method which has been 
used in this work is the truncated Newton method in 
the TNBC package developed by S. G. Nash 4 . The 
software is designed for minimization with bound 
constraints and can be obtained from netlib. Other 
data is necessary to position the airfoils in the proper 
location in three dimensional space. From the airfoils 
in Figure 2, the spanwise locations, h the locations 
of the trailing edges, x e and y e , and the twists, 0 , 
can be determined. It is assumed that each airfoil 
along the fuselage and wing sections are contained in 
a 2 =constant plane and that both airfoils defining 
the winglet are in a y =constant plane. Once the 
airfoils are defined in space, the wing surfaces can 
be constructed to span the airfoils and create the 
airplane surface. 

3. WING SURFACES 

The RAPID methodology creates wing surfaces in 
xyz- space in the form 

X = (x(£, ri), y(£, r]), z{ f, rj)), 0 < £ < 1, 0 < 77 < 1 

The coordinates of surface points are obtained by 
solving the fourth order partial differential equation 


example, if the airfoil generated in Section 2 corre- 
sponded to £ = 0, then the Dirichlet boundary con- 
dition would be 

( x e + s(0 cos( 6 ») + 27 (C) sin( 6 >) 

D 0 (0 = < x e - x(t) sin (0) + y(£) cos{0) 

( h 

The Neumann boundary conditions influence the 
shape of the wing surface and must be determined 
in part from the CAD data. From a geometric point 
of view, the Neumann boundary conditions can be in- 
terpreted as a tangent ribbon which defines the slope 
of the wing surface at the airfoil section. This tan- 
gent ribbon, as depicted in Figure 3, is constructed 
by offsetting the cross sections in Figure 2. For each 
airfoil, an offset distance in the parametric variable 77 
is selected. A new airfoil is constructed at the offset 
location using the RAPID parameters which will give 
a good approximation of the actual CAD cross sec- 
tion at the same location. Although there are a large 
number of parameters which determine the shape and 
location of an airfoil, for calculating the offset airfoil 
it is sufficient to consider only changes in the trail- 
ing edge, (x e ,y e ), chord length, (7, and thickness, T. 
Changes in other airfoil parameters are neglected in 
defining the tangent ribbon. The offset airfoils now 
determine the changes in x and y with respect to 
2 , that is, dx/dz and dy/dz at each point. Now if 
dz/di 7 is given, then dxjdr\ and dy/dy can be cal- 
culated from dx e /dz , dy e /dz , dC fdz, and dT/dz by 
constructing the offset airfoil and applying the chain 
rule. Thus, the procedure for generating the Neu- 
mann boundary conditions can be described in the 
following steps: ( 1 ) calculate airfoil parameters for 
the offset airfoils with chord, thickness, and trailing 
edge coordinates given by 

dC dz dT dz 

+ dz dr] J + dz dr] J 


a 


2 


on 2 

d dr] 2 


X = 0, 


(2) 


where a is a constant weighting parameter. Dirichlet 
and Neumann boundary conditions are imposed at 
77 = 0 and 77 = 1. The boundary conditions are: 


dx e dz 

x e + 

OZ Of] 


Ve + 


dy e dz 
dz dy ’ 


resectively; ( 2 ) construct the offset airfoil with coordi- 
nates (x 0 (£), y 0 (£)); and (3) set the Neumann bound- 
ary conditions 


X(£, 0) = D 0 (£), X(£, 1) = Di(£), 

X„(£,0) = N o (£), X,(^1) = N 1 (0. 

The solution is periodic in £. 

The Dirichlet boundary conditions, defined by the 
functions D 0 (£) and Di(£), are determined from the 
airfoils which form the edges of the surface. For 


( x o (0 — x(£, 0 ) 

N 0 (0 = < Vo(0-V({, 0) 

dz 
V dr] 

The spanwise changes in trailing edge coordinates, 
chord length and thickness can be easily estimated 
from the CAD data. However, an estimate of dz/dr], 
which gives the spanwise influence of the Neumann 
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boundary condition, is not apparent from the data. 
Thus some trial and error was necessary to obtain 
the best fit to the CAD surface. In most cases, 
setting dz/dy equal to the span of the wing section 
is a good starting value. 

Once the boundary conditions are determined, 
equation (2) is solved by generating a Fourier series 
expansion. An approximation of the CAD surface 
using the RAPID technology appears in Figure 4. 

The total number of parameters will now be sum- 
marized. There are a total of six airfoils defining the 
BWB geometry. Each airfoil is defined by eight pa- 
rameters 

C, F, Fi, F 2 , Mi, M 2 , Pi, and P 2 . 

The location and orientation of each airfoil in space 
must be defined. This is done by specifying the span- 
wise location of the airfoil, /i, the coordinates of the 
trailing edge, x e and y e , and the twist angle, 0. There 
are additional parameters for each airfoil used to de- 
fine the Neumann boundary conditions. These are 

dz dx e dy e dC dT 

1 ~ o 1 "~n 5 5 ana. ~ . 

orj oz oz oz oz 

Thus, there are a total of 17 parameters for each 
of six airfoils. Now deleting the spanwise location 
of the centerbody airfoil which we assume to be at 
z = 0, the BWB airplane, as defined by RAPID, is 
generated with a total of 101 parameters. 

Since considerable effort has gone into attempting 
to approximate a given CAD surface with a para- 
metric model, some comparisons of the original CAD 
surface and the RAPID surface will be presented. 
It is evident from comparing Figures 1 and 4 that 
the gross geometry of the BWB airplane has been 
captured in the RAPID model. Figures 5 and 6 
compare the models at two cross sections used in 
defining the RAPID model. The approximation 
of the supercritical airfoil section in Figure 6 is 
excellent, but there is noticeably more error in the 
approximation of the centerbody airfoil in Figure 5. 
A comparison of the planforms appears in Figure 7. 
Whether this degree of fidelity in the approximation 
is sufficient would of course be a judgement call of 
the designer. Further accuracy could be achieved 
with the RAPID technology, but at the cost of 
increasing the complexity of the model and the 
number of parameters. 


4. CFD ANALYSIS 

The RAPID system can output surface grids with 
specified grid point distributions in each coordinate 
direction. These surface grids can be used directly in 
a potential flow analysis code or they can be used as 
boundary surfaces to generate a volume grid for an 
Euler or Navier-Stokes calculation. Similar volume 
grids were constructed for the region about the CAD 
and RAPID models. In order to guarantee com- 
parable grids, a grid was first generated about the 
RAPID surface and then this grid was distorted to 
generate the grid about the CAD surface. The grid 
distortion process was carried out as a application 
of the Coordinate and Sensitivity Multidisciplinary 
Design and Optimization (CSCMDO) system 5 . 
With this system, the differences in the volume 
grids is no greater than the differences in the surface 
grids. The CSCMDO system could also be used 
to generate new volume grids after making small 
changes in design parameters. Solutions to the Euler 
equations for compressible flow were computed on 
grids for the RAPID and CAD geometries. The 
following comparisons are for a free stream Mach 
number of 0.85 and an angle of attack of 2 degrees. 
Figure 8 contains comparisons of contour plots of 
the local Mach number on the upper surfaces of the 
CAD and RAPID geometries. The similarity in the 
solutions is evident, especially the shock structure 
along the wing segment. Further comparisons reveal 
the effects of the differences in geometry. Figures 9 
and 10 are pressure coefficient plots near the center 
of the fuselage cross section and on the wing. The 
pressure coefficient plots are superimposed over 
plot of the airfoil cross sections, so that difference 
in solutions can be compared with the differences 
in geometries. It was interesting to note that the 
RAPID surface geometry gave smoother plots and 
a larger lift coefficient than the solution computed 
with the CAD geometry. 

5. NURBS APPROXIMATION 

NonUniform Rational B-Splines (NURBS) have 
become a standard surface representation in the 
CAD industry 6 . For this reason, a method for 
accurately representing the RAPID surface geometry 
by a NURBS surface has been investigated. As with 
the RAPID methodology, one of the advantages in 
a NURBS representation of a surface is the fact 
that the surface is defined by a set of parameters 
which can be manipulated to change the shape of 
the surface. However, there is a difference in the 
two representations. The RAPID parameters are 
actual geometric lengths and distances, whereas the 
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NURBS parameters are the control points, weights, 
and knots. 

A NURBS surface of degree p in the indirection 
and degree q in the ^-direction is a piecewise rational 
function of the form 

jy, x J2i=0^2j=0 ^i,p( U )Nj,q(v) w iJPiJ 

P (U,V) = 7 —j 

Si=0 Sj=0 Ni,p(u)Nj^ q (v)Wij 

where P ij are the points of the control net, wij 
are the weights, and Ni^ p (u) and Nj^ q (v) are the B- 
spline basis functions defined on the knot vectors 
^o, ' ' ' , fU+p+i and v 0 , • • • , vj+ q +i with 

u 0 = ••• = u p = 0 = % = --- = v q+ 1 

Ul + 1 = • • • = Ul+p+i = 1 = Vj+I = • • • = Vj+q+i. 

Let X m?n , m = 1, • • • , M, n = 1, • • • , N be a sur- 
face grid generated by RAPID. Surfaces as NURBS 
are defined on the unit square and can be evaluated 
at an array of points ( u m ,v n ), m = 1, • • • , M, n = 
1, • • • , N in the mnplane. The optimal NURBS ap- 
proximation is obtained by selecting the control net, 
weights, and knots which yields the minimum value 
of 

M N 

^ ^ £[X m , n P('U m ,'C n )] . 

m = 1 n= 1 

The optimal NURBS will depend on the selection of 
the array of parametric variables u and v. Both uni- 
form and arclength parameters have worked well in 
practice. 

An example of a NURBS approximation of the 
BWB airplane and the control net is plotted in 
Figure 11. This surface is defined using piecewise 
cubic basis functions in both directions and an 8 x 8 
control net. The NURBS surface is therefore defined 
by 64 values for each coordinate of the control net, 
64 weights, and four knot values for each direction. 
Thus the NURBS surface is defined with a total of 
264 parameters. 

The NURBS approximation requires the solution 
of a nonlinear optimization problem. The solution 
is computed using a two-stage process. First, the 
optimal control net is constructed using unit weights 
and uniform knots. This is a linear least-squares 
problem that is solved directly using a QR matrix 
factorization. The resulting set of control points with 
unit weights and uniform knots are initial values for 
calculation of the optimal NURBS. Next, the full 
nonlinear problem with variable control net, weights, 


and knots is solved using the TNBC optimization 
routine. 

6. SUMMARY AND CONCLUSIONS 

The RAPID methodology has been used to 
generate a parametric model for the BWB subsonic 
transport. With only a small set of parameters, there 
is a limit to the accuracy in the approximation of the 
original CAD model. In examining the differences 
in geometry and CFD solutions with the CAD and 
Rapid models, it should always be noted that this 
methodology is intended for aerodynamic analysis 
at the conceptual stage of the design process. At 
this stage, the CAD geometry may only be a crude 
approximation of the final configuration. The ob- 
jective in constructing the RAPID model is to have 
a parametrization of the BWB configuration. The 
RAPID model can then be used in a parameter study 
to predict the effects of changes in the geometry and 
improve upon the original design. 
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